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H 

£jj ' The 7-ray and neutrino emissions from dark matter (DM) annihilation in galaxy 



clusters are studied. After about one year operation of Fermi-LAT, several nearby 
clusters are reported with stringent upper limits of GeV 7-ray emission. We use the 



CO 
> 

\ Fermi-LAT upper limits of these clusters to constrain the DM model parameters. 



• We find that the DM model distributed with substructures predicted in cold DM 

(CDM) scenario is strongly constrained by Fermi-LAT 7-ray data. Especially for the 
leptonic annihilation scenario which may account for the e ± excesses discovered by 
PAMELA/Fermi-LAT/HESS, the constraint on the minimum mass of substructures 
is of the level 10 2 — 10 3 M , which is much larger than that expected in CDM picture, 
but is consistent with a warm DM scenario. We further investigate the sensitivity of 
neutrino detections of the clusters by IceCube. It is found that neutrino detection is 
much more difficult than 7-rays. Only for very heavy DM (~ 10 TeV) together with a 
considerable branching ratio to line neutrinos the neutrino sensitivity is comparable 
with that of 7-rays. 

PACS numbers: 95.35.+d,95.85.Pw,95.85.Ry,98.65.Cw 
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INTRODUCTION 



The existence of dark matter (DM) has been established by many astrophysical observa- 
tions, but the nature of DM paticle is still unclear. Among the large amount of candidates 
proposed in many theories of new physics, the weakly interacting massive particle (WIMP) 
is the most popular and attractive one [ll, |2[ . The mass of WIMP is generally from a few 
GeV to TeV, and the interaction strength is of the weak scale, which can give the right relic 
density of DM. In this scenario, the weak interaction of DM particles would produce observ- 
able standard model particles, such as charged anti-matter particles, photons and neutrinos. 
Investigating such particles from the cosmic rays (CRs) is the task of DM indirect detection. 

The recently reported new signatures of CR positrons, antiprotons and electrons by 
PAMELA [3], Q, ATIC fl, HESS |s, 7| and Fermi-LAT jg] have stimulated great interests 
and extensive studies of the DM indirect searches. The DM scenario with mass O(TeV), 
leptonic annihilation/decay final states and a high annihilation/decay rate can well explain 
the observational data (e.g., js), [lo|). Furthermore, more quantitative constraints on the DM 
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12|. 



model parameters can be derived through a global fitting method 

Regardless of detailed models of DM to explain the data, it is essential to find observable 
signals to test the models. Since the charged particles will gyrate in the magnetic field and 
lose most of the source information, it is difficult to test the DM models using only the data 
of charged CRs. Gamma-rays and neutrinos seem to be very good probes. There are several 
advantages of using 7-ray photons and neutrinos to investigate the DM models. Firstly, 
photons and neutrinos propagate along straight line and can trace back to the source sites 
where the DM annihilation/decay takes place. Secondly there is little interaction during the 
propagation and most of the primary source information hold. Thirdly the effective volume 
of which photons and neutrinos can reach is much larger than that of charged particles, e.g., 
from the Milky Way to extragalactic space, and even the early Universe. It has been shown 
in some works that 7-rays and neutrinos can be powerful tools to test the DM scenarios 
explaining the CR lepton data (e.g., jl3l-l20|. |2l|-|27l|). 



There are many sites proposed to be goo d candidates for the search o 



from DM, such as the Galactic center 
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Galactic halo 
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20 



7-rays and neutrinos 



or substructures 29M32I] . the extragalactic space 33M35I] and the emissions at the early 



28[, satellite galaxies 



Universe 



3(3 



. As the largest gravitational bounding system in the Universe, galaxy 
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clusters may also be useful for DM indirect searches [39|. Pinzke et al. investigated the 
7-ray emission from nearby clusters and used EGRET upper limits to set constraints on 
the DM model parameters 40|. They found that if the DM annihilation was responsible 
for the electron/positron excesses and the luminosity- mass distribution of DM substructures 
in clusters followed the extrapolation of numerical simulation results, the minimum mass 
of DM subhalos should be larger than 10 -2 M Q in order not to exceed the EGRET limits. 
This is a useful way to study the particle nature of DM through structures. 

After more than one year's operation, Fermi-LAT reported some results about the 7-ray 



emission from galaxy clusters 



411 ] . Non detection of significant 7-ray emission from galaxy 



clusters was reported except for Perseus cluster, in which the emission from the central 
galaxy NGC 1257 was discovered {42]. The upper limits given by Fermi-LAT are lower by 
more than one order of magnitude than that given by EGRET. It can be expected that 
the new results from Fermi-LAT will set much stronger constraints on the DM models. 



In Ref. 



4ll | the constraints on DM mass and annihilation cross section were presented 



assuming fi + fi and bb channels. In this work we will also use the Fermi-LAT upper limits 



4l| . we will pay more attention 



to constrain DM model parameters. Different from Ref. 
on the implication of DM structure properties such as the minimal mass of subhalo M m i n , 
which would be important for understanding the nature of DM particle. This is one of the 
motivations of this study. 

Another motivation of this work is the neutrino emission. Neutrinos can be served as 
an independent diagnostic of DM indirect searches besides photons. It has been shown 
that the measured atmospheric neutrino background can set effective constraints on the DM 



annihilation cross section 



43 



441 ]. There are no high energy astrophysical neutrinos being 



detected currently, so it is valuable to explore the sensitivity of the forthcoming neutrino 
detectors to the neutrino signals from the DM annihilation. Due to the very weak interaction 
cross section between neutrinos and matter, we generally need large detector volume. The 
ongoing experiment IceCube has an effective volume ~km 3 , which would give unprecedented 
sensitivity for the neutrino detection up to very high energies. The detectability of neutrino 
signals from DM annihilation in galaxy clusters by e.g. IceCube, will be discussed in this 
work. 

This paper is organized as follows. In Sec. II, we discuss the 7-ray emission from several 
galaxy clusters, and employ the recent Fermi-LAT limits of these clusers to constrain the DM 



i 



model parameters. In Sec. Ill, we discuss the detectability of neutrino emission from the 
galaxy clusters by the neutrino detectors. The last section is our conclusions and discussions. 

II. GAMMA RAYS FROM GALAXY CLUSTERS 
A. Cluster sample 

It is known that the objects with high masses and small distances will be very efficient for 
the DM searches. Therefore nearby massive clusters are the first choice of study. Here we 
adopt a sample of 6 clusters with redshift from 0.0031 to 0.0231 (corresponding to distance 
from 13 to 100 Mpc for a standard ACDM cosmology), which are reported with flux upper 
limits by Fermi-LAT. The basic parameters of these clusters are compiled in Table [B 

The flux of DM annihilation from cluster is generally scaled with M^o/d 2 , where a 
depends on the concentration-mass relation and DM profile of the halo. In this work we will 
assume NFW profile for the halo of cluster. The final results are expected not sensitively 
dependent on the profile since most of the cluster lies in the angular window of Fermi-LAT 



(3.5° for 100 MeV, 45j). As a benchmark configuration, we adopt the concentration-mass 



relation fitted from X-ray observations 



46) 



9.0 / M vW 



1 + z \10 u h- 1 M Qi 
After correcting the definition of virial overdensity in Ref. 

have [47] 



(1) 



46| (A « 100) to A = 200 we 



6.9 ( M 200 y°- 178 

C200 " TT~z x ^lo^ ) ■ (2) 

For this concentration-mass relation we find a ~ 0.65. Comparing the quantity Mf 00 /cP 
among these clusters, we find that DM signals from clusters NGC 4636, M49 and Fornax 
are of the same level, and are several times larger than the rest three clusters. In the 
following we will see that these three clusters will indeed give stronger constraints on the 
DM models. 

B. Gamma-ray emission from DM distribution in clusters 



There are generally two kinds of 7-ray emission from DM annihilation: one is produced 
directly from the annihilation final state particles which is called primary emission (such as 
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TABLE I: Cluster sample 



Name 


z a R.A." 


Dec. Q 


m 2OO (io 14 m ; 


6 r 200 (Mp C ) 6 


NGC 4636 0.0031 12 h 43 m 


2°41' 


0.25 


0.60 


M49 


0.0033 i2 h m m 


8°00' 


0.46 


0.73 


Fornax 


0.0046 03 h 39 m 


-35°27' 


1.00 


0.95 


Centaurus 


0.0114 12 h 49 m 


-41°18' 


2.66 


1.32 


AWM 7 


0.0172 02 h 55 m 


41°35' 


4.28 


1.54 


Coma 


0.0231 13 h 00 m 


27° 59' 


13.65 


2.27 



a Redshift and coordinates are adopted from NASA/IPAC Extragalactic Database, 



|http://nedwww.ipac. caltech.edu/ 



b Virial mass and radius parameters are taken from Ref. 



the 7 rays by 7r° decay after hadronization or emission directly from final charged leptons), 
and the other is produced through interactions of final state particles with external medium 
or radiation field such as the inverse Compton (IC) radiation which is called secondary 
emission hereafter. The primary 7-ray flux observed on the Earth from DM annihilation in 
a galaxy cluster can be expressed as 

_ M dN fp\r)dV 

^ ~ 2ml dE 4ird? L ' [6) 

x >-> 

where m x is the mass of DM particle, (av) is the annihilation cross section of DM, ^ is the 
yield spectrum of 7-rays per annihilation which is simulated using PYTHIA j^], di is the 
luminosity distance of the cluster, p(r) is the density distribution of DM inside the cluster 
with r the distance from the cluster center. All of the cluster is taken into account in the 
integral since the analysis of Fermi-LAT was done in a 10 degree radius of each cluster 41], 



which is large enough to contain the whole cluster halo. For the smooth halo we assume the 
density distribution to be NFW profile [so | 

Psm( r ) = 7 / ryz : 1 (4) 
(r/r s )(l + r/r s y 

where parameters r s and p s can be determined by the concentration-mass relation and 
normalization of total mass. 

Since there are substructures in the clusters, such as galaxy groups and galaxies, we have 
to take these into account. We will see later that the existence of substructures enhances 
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the annihilation luminosity of DM and is the main reason that affects the 7 flux. To take 
into account the effect of substructures, we replace p 2 in Eq. (J3J) with p 2 ot = p 2 m + (p 2 uh ), 
where the average density square of substructures reads 

(pUr)) = J *M^m X L(M) ' (5) 
in which d y^ M is the number density of subhalos in mass bin dM, L(M) = J v p 2 uh dV is 



the intrinsic annihilation luminosity of a subhalo with mass M. 



n this work we will employ 



the results from recent high resolution simulation, Aquarius [Ml, |52[ to treat the subhalos. 
Because the concentration and density profile of subhalos are very complicated inside the 
host halo, the detailed computation using Eq. (j3J) is difficult. Thus we directly adopt the 



counted results of luminosity distribution from the simulation 1 53] 



d£ ^ , _ x -0.1 , _ x -2.9 



r M) = x L(M) oc 1 + x M b . (6) 

dikT ' ; d^dM 1 ; V0.27W V 0.2r 200 ; 1 ; 

Similar with Ref. 40[ we adopt a scale between the Milky Way like halo given in Aquarius 
simulation and the case of clusters, i.e., the ratio of £ su t,/£ sm keeps unchanged whatever the 
mass is. The maximum mass of subhalos found in simulation is about 0.01Mh ost . But the 
minimum mass is not well known due to the limit of resolution of the numerical simulation. 
From the observational point of view, we have observed DM halos with mass ~ 10 7 M , e.g. 
dwarf galaxies. While the study of free streaming of cold DM (CDM) particles indicates 



a minimum halo mass down to ~ 10 7 M m 



541 ] . In this work we leave M m ; n to be a free 



parameter and investigate the effects of M min on the DM signals. 

Besides the primary 7-ray emission, there is also secondary production of 7-ray pho- 
tons through the IC scatterings between the DM induced electrons/positrons and the cos- 
mic microwave background (CMB) field. For the calculation of the fluxes of secondary 



IC emission please see the Ref. [I6j|. Note that when calculating the energy loss rate of 
electrons/positrons, both the IC loss induced by scattering with CMB photons and the syn- 
chrotron loss in the magnetic field in clusters are considered. The average value of magnetic 



field strength is assumed to be ~ 1 pG [55]. 



1 This relation is different from that given in Ref. [5l| . where L(> M) m M 226 was found. According 
to this fit we have L(> M) rs (M~°' le — M~^ 6 ), which gives similar behavior as that in Ref. 51 



the large (resolved) mass range of the simulation, but is different when extrapolating to low (unresolved) 
mass range. 
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FIG. 1: Integral spectra of the IC and FSR components from DM annihilation in Fornax cluster. 
The left panel is for channel, and the right panel is for bb channel respectively. Also shown 

are the 95% upper limits (arrows) of Fermi-LAT 1 yr observations. See the text for details. 



For illustration we show the integral spectra of 7-rays from DM annihilation in Fornax 
cluster in Fig. [TJ Here we adopt a sample DM model with m x = 1 TeV, (crv) = 10 -23 
cm 3 s _1 , and the annihilation channels with (left) and bb (right) respectively. The 

primary and secondary components are shown separately. In each group we show three 
curves which represent the smooth halo contribution, the total emission with subhalos down 
to two different M min . The 95% confidence level upper limits from Fermi-LAT are shown 
by arrows. It is shown that in the energy range interested here, i.e. 0.1 — 10 GeV, the 
secondary radiation from IC is dominant for fi + fi~ channel, while for bb channel the primary 
contribution is dominant. This is because for final state the spectra of photons and 

electrons are relatively hard, and the secondary produced photons through IC can just lie 
in the interested energy range. For bb final state the energies of photons and electrons from 
the hadronic cascade are generally much lower than the mass of DM, so the IC component 
dominates at even lower energies (< 0.1 GeV). This conclusion will always hold for m x 
ranging from 100 GeV to 10 TeV for fi + fi~ channel (10 GeV to 1 TeV for bb), which is the 
sensitive region explored by Fermi-LAT. 

In 41] the Fermi-LAT has derived preliminary constraints on the 7-ray fluxes from DM 
annihilation from these clusters by assuming the and bb final states. According to 

these results, we derive the constraints on the m x — (crv) plane. In Figs. [2] we show the 
constraints for the /x 4 ' \i~ channel. It is shown that the constrains of the smooth component 



s 



NGC4636 l0 9 10 M min M49 lo 9io M , 




FIG. 2: Fermi-LAT constraints on the DM model parameters m x and (av) for different minimum 
halo mass M m j n . The DM annihilation channel is ■ Dashed circles are the 3a and 5a param- 

eters regions which can fit the PAMELA/Fermi-LAT/HESS data of the CR positrons/electrons 
56]. 

(that M min = 1O 12 M ) is still very weak. If the substructures are taken into account the 
constraints can be stronger by more than one order of magnitude, depending on the free 
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FIG. 3: Comparison of the constraints 
Via Lactea substructure model (right, 



:'or Maccio et al. (2008) concentration model (left, [571]) and 



streaming mass of DM subhalo. If the minimum mass of DM subhalos can be as low as 
ICG 7 — ICG 6 M as predicted in the supersymmetric DM scenario, the constraint on the cross 
section of TeV DM can reach ~ 10~ 24 cm 3 s -1 . This constraint is actually powerful enough 
to explore the DM models which are proposed to explain the CR lepton excesses. In Fig. 
El we plot the favored parameter regions of the DM model to explain the PAMELA/Fermi- 



LAT/HESS data of the CR positrons/electrons with final state 56j. We can see that 

for almost all of the 6 clusters the PAMELA/Fermi-LAT/HESS favored parameter regions 
can be excluded if M min is down to ~ ICG 6 M Q . If DM annihilation is responsible for 



the observationa 



positron/electron excesses, it may indicate that the cutoff of M min should 



be much larger [40]. For clusters Fornax and NGC 4636, the constraint of M m \ n is about 
10 2 - 10 3 M for the best fit mass and cross section 2 . This constraint is much stronger than 
that derived using EGRET upper limit about Virgo cluster j^ . 

This cutoff of the mass of substructures may have important implication of the particle 
nature of DM. It gives an estimate of the free streaming scale of the matter power spectrum 
as k < 750 Mpc" 1 , which is not very far from the lower limits given by Lyman-a power 



spectrum measurements 6l|, [62|. Compared with the canonical value expected in CDM 



2 Note that here we assume a constant boost factor of the DM models responsible to the electron/positron 
excesses. However, as proposed in the literature the large boost factor might be due to velocity-dependent 



annihilation cross section such as the Sommcrfcld effect 



59 



60] . In that case the boost factor might be 



different in clusters from in the Galaxy. We do not discuss this effect in detail in this study. 
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FIG. 4: The same as Fig. [2] but for bb final state of DM annihilation. 



picture, it favors a warm massive DM scenario which may be produced non-thermally in the 



early Universe 



63l | . Such nonthermally produced DM particles have large initial velocity 



and large free streaming. Thus the matter power spectrum is suppressed at small scales and 
leads to less low mass subhalos k34, 



65|. 



As a check of the model uncertainties of the halo structure configuration, we compare 
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the results using Maccio (2008) concentration-mass relation given in Ref. 57 ] 

-0.098 



c 200 



3.56 
1 + z 



X 



200 



10 15 M m 



(7) 



For substructures we also use the results from another high-resolution simulation Via Lactea. 
According to Fig. 3 of Ref. 58], the substructure enhancement is simply extracted to be 
B = 10 x (M"^ 048 - M"° x 048 ) with M max = 0.01M host . Note that there is a host halo mass 
dependence of B as given in Ref. 58j. However, since the results with Mh ost larger than the 
galaxy scale halo were not calibrated in the simulation, we adopt B{M\ iOSt ) at Mh os t = 10 12 
M Q and apply it to the cluster scale halos. This treatment is also consistent with the 
above assumption that we adopt the scaling relation of Eq. to be the same for Milky 



Way halo and cluster ha 
dependence of B in Ref. 



os. The extrapolation to the cluster masses following the Mh os t 



5S 



would lead to a two times larger boost factor. The constraints 
for channel from Fornax cluster for these two models are shown in Fig. [3j It is 

shown that for Maccio (2008) concentration the constraint is about 2 times weaker than the 
benchmark model. And for Via Lactea substructures the total substructure enhancement 
is some weaker than that of Aquarius simulation. However, in both of these cases we see 
that the PAMELA/Fermi-L AT/HESS favored parameter regions can be constrained, given 
M m i n ~ 10 -6 M . For the best fit mass and cross section the constraint on M min is about 
10° - 10 2 M . 

The results for bb channel are shown in Fig. HI We can see that for m x ^ 100 GeV 
the most stringent constraint of (av ) can reach the thermal scale of 3 x 10~ 26 cm 3 s" 1 for 
M m i n = 10 -6 M Q . The future observation of Fermi-LAT can put a strong constraint on the 
DM model even for the thermal scenario. Note that for bb channel only the primary emission 
is assumed when deriving the flux limits. If the IC component is included the constraints 
will be stronger for relatively heavy DM mass (~TeV). Finally we should point out that bb 
channel is typically not suitable to explain the lepton excesses observed by PAMELA/Fermi- 
LAT/HESS, due to both the constraint from PAMELA p/p data j^, 9] and the spectral shape 
required by lepton excesses 10|. Here we include the study of bb channel is just to show the 
power of Fermi-LAT to the supersymmetric-like DM particles. 
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III. NEUTRINOS 

In this section, we discuss the sensitivity of detecting neutrino signals from clusters. 
The cluster could be treated as high energy neutrino point source, and it is possible to be 
observed at the on-going large volume neutrino telescopes, such as IceCube. To suppress 
the large atmospheric muon background, the neutrino telescopes usually detect the upward 
muons induced by muon neutrinos through interacting with the matter surrounding the 
detectors. Therefore, the south pole based detector IceCube is more suitable to probe the 
neutrino sources in the northern hemisphere. For the sample of clusters in Table [J Fornax 
and Centaurus locating in the southern hemisphere are not good candidates of Icecube 3 . 
AWM 7 with declination of 41°35' and Coma with declination of 27°59' are suitable for 
IceCube, but such two clusters are more distant away from us than other clusters. Taking 
into account the location, mass and distance, we find M49 with declination of 8°00' is better 
to be detected than others. 

Similar to the primary 7-ray flux, the neutrino flux from DM annihilation in the cluster 
can also be calculated according to Eq. fl3]), by replacing ^| with ^\ u - In the following 
we will mainly discuss two DM annihilation channels, fi + fi~ and fi + fi~ + v^p^. The bb 
channel as discussed in the previous section will also be mentioned. However, as we will see 
below, it gives negligible neutrino signals. We use the PYTHIA [49] to simulate the initial 
neutrino spectra from decay of annihilation final states. We further assume the neutrino 
flavor distribution is 1 : 1 : 1 at the Earth due to vacuum oscillation during the propagation. 

The through-going upward muon rate at the detector can be calculated as 



dN„ f ,„ r x ,„ dN, 



dQ / dE„ 



dE^ J J E v>1 dE v 



R(E.) + (u^D), (8) 



- n p + [p — > n 



where n p (n n ) is the number density of protons (neutrons) in the matter around the detector, 
R{Efj) named muon range is the distance that a muon can travel in matter before its energy 
drops below the threshold energy of detector E t h_, which is given by 



3 The IceCube + DeepCore has the capability to search the downward neutrinos, but the angular resolution 
is fairly bad j(36j]. It is very difficult to distinguish the high energy neutrino sources from the high 
atmospheric neutrino background without powerful angular resolution. 
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with a, /3 the parameters describing the energy loss of muons as dE^/dx = —a — jSE^. 

The main background for high energy neutrino detection is the atmospheric neutrinos. 
The atmospheric neutrino flux decreases rapidly as energy increasing. We use a parametriza- 



tion of atmospheric neutrino flux 



671 ] which describes the results of Ref. 



(is 



as 



^ 0.018 , 0039 



dE u dtt v \l + 0m4E„\cos9\ 1 + 0.00139^| cos^| 

where N = 1.95 x 10 17 (1.35 x 10 17 ) GeV _1 km _2 yr _1 sr _1 for v,j, (z/ M ) respectively, 9 is the 
zenith angle. 

IceCube could have an angular resolution ~ 1°, which is effective to suppress the diffuse 
atmospheric neutrino background. Notice the cluster is not an exact point source, we uti- 
lize an angular resolution of 3.0° (1.5°) for M49 (AWM47/Coma) cluster. The number of 
atmospheric neutrinos for 3° resolution angle is ~ 4 times larger than it for 1.5° cone. 
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FIG. 5: Muon flux induced by neutrinos from DM annihilation in M49 cluster, for muon (left) 
and bb (right) final states respectively. The mass of DM is 1 TeV and cross section is 10 -23 cm 3 
s _1 . The results of smooth halo (green) and smooth halo together with subhalos with two values 
of M m j n , 10 -1 M (blue) and 10 -6 M (red) are shown. In the left panel, the solid lines denote 
channel, and dashed lines denote plus Vupn with equal branching ratios. Also shown 

is muon flux induced by atmospheric muon neutrino background in 3° cone. 



In Fig. we show the through-going muon flux induced by DM annihilation in M49 
cluster. Similar as in Fig. [T]we adopt m x = 1 TeV, (av) = 10~ 23 cm 3 s _1 , and annihilation 
channels are (left solid), + v ll y >1 (B^ = B v = 0.5, left dashed) and bb (right). For 
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bb channel, the neutrinos are produced through decay of hadrons induced by bb hadroniza- 
tion. It is shown that the muon spectrum of such channel is very soft and difficult to detect 
given the high level of background. The case for channel seems better but the signal 

is still very weak. Even for M m i n = 10~ 6 M , the muon flux from DM annihilation is ~ 100 
times smaller than the background in energy range (200, 800) GeV. Only for fi + fi~ + v^v^ 
the situation is better. The reason is that monochromatic neutrino spectrum is harder than 
other channels, and is easier to be detected. Compared with Fig. (TJ it is not strange to see 
that the neutrino detection sensitivity would be much weaker than the 7-ray detection. 
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FIG. 6: (av) required to discover neutrinos from DM annihilation in M49 cluster as a function 
of DM mass. The left panel is for channel, and the right panel is for /U + ^/~ + v^v^ with 

branching ratios = B v = 0.5. Dashed circles are the 3a and 5a parameters regions which can fit 



the PAMELA/Fermi-LAT/HESS data of the CR positrons/electrons [56(]. The dot-dashed curves 
represent the 7-ray constraints for M m - m = 10 _1 and 10 -6 Mq respectively (see Fig. [2]). The circles 
and 7-ray constraints in the right panel are scaled upwards by a factor 2 due to the branching ratio 



B„ = 0.5. 



The total muon event rate in a specific energy bin at the detector is 

dA^ 



N 



dE^A cS (E„6)AT, 



(11) 



where A e g is the effective muon detecting area taken from Ref. 69|, AT is the operation 



time which is set as 10 years here. We take the threshold energy to be E t ^ ~40 GeV, and 



assume the energy resolution is Alog 10 -E = 0.35 



66|. 
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FIG. 7: The same as Fig. [6] but for AWM7 cluster. 
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FIG. 8: The same as Fig. [6] but for Coma cluster. 



In Fig. [6] we give the IceCube sensitivity of detecting neutrinos from DM annihilations 
in M49 for 10-yr exposure. The left and right panels show the results for annihilation 
channels fi + n~ and + respectively. For channel and M min = 10 -6 M Q , a 5<r 

detection requires (av) as large as O(10~ 22 ) cm 3 s~ 1 . We can see such a cross section is much 
larger than that required to explain PAMELA/Fermi-LAT/HESS data. For the channel to 
equal /i + n~ + v^v^ the sensitivity is ~ 10 times better than channel. Because the 

neutrino- nucleon cross section and muon range are approximately proportional to E v , while 
the atmospheric neutrino background decreases as ~ E~ 3 , the neutrino telescope is more 
powerful to detect the high energy neutrinos. For the same reason, the detector is more 
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sensitive to explore heavy DM. It is shown that if DM annihilation products have a large 
branching ratio to vv, IceCube could reach the parameter space to explain PAMELA/Fermi- 
L AT/HESS results. 

The sensitivities for other two clusters, AWM 7 and Coma, are shown in Figs. [7] and [8] 
respectively. These two clusters are much more distant from us than M49. Note for these 
two clusters we take 1.5° cone, which is enough to include all of the cluster, to count the 
atmospheric background. Due to a much lower level of background and larger masses of 
AWM 7 and Coma, the sensitivities are only several tens percent weaker than M49. Because 
the declinations of AWM 7 and Coma are larger than M49 which is close to the horizon, 
it would be more effective to reject the CR muon background and could provide clearer 
detection of signals. 

Compared with the 7-ray sensitivity discussed in Sec. II, the sensitivity of neutrino 
detection is relatively poor. If we employ the 7-ray constraint of DM annihilation in clusters 
(e.g., M min » 10 3 M for PAMELA/Fermi-LAT/HESS best fit parameters), there would be 
little chance to detect the neutrino signals. Only for very heavy DM (m x ~ 10 TeV), the 
sensitivities of the two ways are comparable. 

IV. SUMMARY AND DISCUSSION 

In this work we investigate the 7-ray and neutrino emission from DM annihilation in 
clusters of galaxies. A sample of several nearby clusters is considered. Both the annihilations 
in the host halo and substructures are taken into account. For the annihilation luminosity of 
substructures we adopt the result of Aquarius simulation, and scale it from Milky Way like 
halo to the cluster like halo. There is also a component of unresolved substructures which is 
not seen due to the limit of resolution of numerical simulation. For the contribution from the 
unresolved substructures we adopt an extrapolation of the luminosity-mass relation fitted 
from the resolved subhalos in the simulation. The minimum mass of the subhalo M m i n is left 
to be a free parameter. The value of M m ; n may catch the information about free streaming 
length, and may reflect the generation history of the DM particle. 

For the 7-ray emission we consider two typical annihilation channels, quark final states 
bb and lepton final states Both the primary component of photons from hadron 

decay and final state radiation, and the secondary component which is produced by the IC 
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scattering of DM-induced electrons/positrons and CMB field are taken into account. The 
Fermi-LAT upper limits of the 7-ray emission from these clusters are employed to constrain 
the m x — (av) parameter plane. The constraints depend on the value of M min . For M min 
down to ~ 10 -6 M , which is expected for the structure formation of neutralino-like CDM 
picture, the constraints are ~ 50 times stronger than the case of smooth halo only. Typically 
for m x = 1 TeV and M min = 10~ 6 M Q the strongest constraint on cross section from the 
cluster sample is 10 -24 (10 -24 ) cm 3 s _1 for (bb). 

It is of great interest for the fi + fi~ channel which is proposed to explain the positron 



and electron excesses reported by PAMELA, Fermi-LAT and HESS [56|, [70] . A very large 
annihilation cross section (or boost factor) is needed to explain the data. It is shown that the 
Fermi-LAT observations about 7-rays from galaxy clusters can strongly constrain the model 
parameters. If we fix the mass and cross section to the values explaining the e ± excesses, 
the minimum mass of substructures M m [ n is constrained to be larger than 10 2 — 10 3 M Q . 
Such a large value of halo mass means a very large free streaming length of DM particle. It 
may indicate the nature of DM particles is warm instead of cold 
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6;: 



Finally we calculate the sensitivity of detecting neutrinos from these clusters by the 
IceCube detector. It is shown to detect neutrinos would be much more difficult than 7-rays. 
For bb final state the sensitivity is extremely poor due to the neutrino spectrum from bb 
hadronization is soft and suffers from a very high atmospheric background. The case becomes 
better for fi + fi~ final state. However, the signal is still very weak. For example, even for 
M min = 10~ 6 M , a 5a detection with ~ 10-yr exposure requires (av) as large as O(10~ 22 ) 
cm 3 s -1 . This result does not have enough capability to explore the PAMELA/Fermi- 
LAT/HESS favored parameter region. If we consider the model with a large fraction of 
line neutrino emission the detectability would be much better (e.g. improved by an order 
of magnitude). However, compared with 7-rays the constraint from neutrinos is still some 
weaker. Only for very heavy DM (m x ~ 10 TeV), the sensitivity of neutrino detection can 
be comparable with 7-rays. 
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